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Abstract — The essential insight of quantum error correction 
was that quantum information can be protected by suitably en- 
coding this quantum information across multiple independently 
erred quantum systems. Recently it was realized that, since 
the most general method for encoding quantum information 
is to encode it into a subsystem, there exists a novel form 
of quantum error correction beyond the traditional quantum 
error correcting subspace codes. These new quantum error 
correcting subsystem codes differ from subspace codes in that 
their quantum correcting routines can be considerably simpler 
than related subspace codes. Here we present a class of quantum 
error correcting subsystem codes constructed from two classical 
linear codes. These codes are the subsystem versions of the 
quantum error correcting subspace codes which are general- 
izations of Shor's original quantum error correcting subspace 
codes. For every Shor-type code, the codes we present give a 
considerable savings in the number of stabilizer measurements 
needed in their error recovery routines. 

Real quantum systems are open quantum systems which 
can couple in an unwanted manner to an environment or con- 
trol system and lose their intrinsic quantum nature through 
the processes of decoherence, quantum noise, and impre- 
cise measurement, preparation, and control. These problems 
serve as an obstacle towards the eventual construction of a 
robust large scale quantum computer[l], [2], [3], [4]. If left 
unchecked, these problems turn a quantum computer into a 
classical information processing device, or even worse, into a 
machine which can enact no computation at all. Fortunately, 
however, soon after the discovery that quantum computers 
were more efficient at solving certain problems than classical 
computers [5], [6], it was discovered that, under specific rea- 
sonable physical assumptions, a fault-tolerant quantum com- 
puter could be built. In particular, a set of threshold theorems 
for fault-tolerant quantum computation were established[7], 
[8], [9], [10], [11], [12]. These theorems proved (or gave 
a heuristic proof) that if decoherence, quantum noise, and 
lack of control were all small enough in comparison to the 
ability to control the quantum system (below some threshold 
or thresholds), then these noisy imprecise bare devices could 
be efficiently put together in a fashion which decreased the 
failure probability of a quantum computer to any desired 
level. 

A central insight used in the theory of fault-tolerant quan- 
tum computation is that quantum information can be encoded 
into what is known as a quantum error correcting code[13]. 
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[14]. Such codes spread quantum information across many 
physical quantum bits (qubits) of a system and protect the 
quantum information so encoded from the undesired effects 
which cause a loss of quantum coherence. In the original 
theory of quantum error correction, quantum information 
was encoded into a subspace of the Hilbert space of many 
quantum systems. Such subspace codings, however, are not 
the most general way to encode quantum information into 
a quantum system. The most general method to encode 
quantum information is to encode it into a subsystem[15]. 
This realization has recently led to the discovery of a 
new method for perform quantum error correction using 
quantum error correcting subsystems [16], [17], [18], [19]. 
While every quantum error correcting subsystem can be 
turned into a subspace code, subsystems codes differ sig- 
nificantly in how quantum error correction is performed on 
the encoded quantum information. Of particular significance 
is that quantum error correcting subsystems can significantly 
reduce the number of stabilizer measurements needed during 
their quantum error recovery routines [20], [21], [22]. This in 
turn can lead to significantly improved thresholds for fault- 
tolerant quantum computation[23], [24]. 

In this paper we construct a class of new quantum error 
correcting subsystem codes. These codes are a generalization 
of a code presented by one of us in [21] where they were 
put to use in an attempt to construct self-correcting quantum 
memories. The codes we describe are in the class of stabilizer 
subsystem codes described by Poulin[20]. In particular they 
can be described as the subsystem version of the class 
of codes arising from generalizing Shor's original quantum 
error correcting subspace codes [13]. The subsystems codes 
we describe are constructed by taking two classical linear 
codes and, instead of concatenating them in the manner 
of Shor, constructing a new quantum code which is a 
subsystem quantum error correcting code directly from these 
classical linear codes. In Shor's original construction a bit 
flip error correcting code is concatenated with a phase flip 
error correcting code (or vice versa.) The order of these 
concatenations presents an asymmetry in the recovery routine 
for these codes. If, for example a bit flip code is used on 
the lowest level, then the recovery procedure for the bit flip 
code must be enacted for every lowest level code, whereas 
the phase flip code recovery routine need only be enacted 



on the next level of the code. In the codes we construct, 
this asymmetry between bit flips and phase flips in the 
recovery routine is removed. This leads to codes which have 
substantially simpler error recovery routines, but provide 
the same amount of protection as the generalized Shor- 
type codes. Since the threshold for fault-tolerant quantum 
computation depends significantly on the complexity of the 
error recovery routine, these codes offer the hope of a 
substantial improvement in the threshold for fault-tolerant 
quantum computation[23], [24]. 

I. Classical and Quantum Linear Codes 

We begin by briefly reviewing classical and quantum 
linear codes in order to estabUsh conventions for our future 
constructions. 

We will work with classical codes over the finite field 
F2. F2 is the field with two elements and 1 with addition 
defined as 0-^0=1-1-1 = and 0-1-1 = 1-^0=1 and 
multiplication = 01 = l- = and 11 = 1. Ordered 
n- tuples of F2 form a Unear space, (F2)". A classical binary 
linear code on n bits is a subspace C of (F2)". A subspace 
can be described by a set of k basis vectors. We can list 
such a set of k basis vectors in a fc x n matrix, G, with 
elements from F2. This matrix is called the generator matrix 
for the code C. The generator matrix defines the encoding 
procedure for the code. In particular if we left multiply the 
kxn matrix G by a fc-tuple row vector of bits to be encoded, 
then we obtain the encoded length n string for these bits. 
For a code C, the parity check matrix, P is an (n — fc) x n 
matrix with elements from F2 such that Pv = for all 
w e C. Note that we use a slightly non-conventional row 
and column ordering for the parity and generator matrices. 
The (Hamming) distance between two elements of F2, v and 
w, is given by H{v,w) = Y^"^i Svi,wi, i-C- the number of 
places where the two n bit strings disagree. We say that a 
classical linear code is a [n,k,d\ code if it encodes k bits 
into length n bits and the minimum distance between any 
two elements of the code space is d. 

We now discuss quantum linear codes which are also 
known as stabilizer codes. We will assume that the reader is 
familiar with stabilizer codes at the level of [25] and provide 
the following review to set our definitions. 

Suppose we have n qubits. Then V is the group, com- 
monly known as the Pauh group, whose elements are 
made up of tensor products of single qubit Pauli operators 
(/, X, Y, Z) along with a global phase of i", for s G 
{0, 1, 2, 3}. Every Pauh operator can be written in the form 
^s_^ai j^bi g, z'^^X''^ (g) • • • (g) Z''"X''" where s G {0, 1, 2, 3} 
and ai,bi G {0,1}. The weight of a Pauli group element 
is the number of non-identity single qubit operators in the 
above expansion. All Pauh group elements either connmute 
with each other [P, Q] = PQ — QP = or anticommute 
wifli each other {P, Q} = PQ + QP = 0. 

A stabilizer group is an abelian subgroup of the Pauh 
group which does not contain the Pauli group element — /®". 
All of the elements of a stabihzer square to /®" and hence 
have eigenvalues ±1. Given a stabiUzer group <S, we can 



define a stabilizer code as the subspace of the Hilbert space 
on n qubits which is stabilized by all of elements in <S, 
i.e. the subspace Hs defined as all states \ip) such that 
Slip) = \tp) for all S G S. If the stabilizer group has a 
minimal set of generators Si, . . . , Sn-k, then the dimension 
of the stabilized subspace (the code subspace) is 2'^, i.e. it 
encodes k qubits of quantum information. The normalizer 
Af of the StabiUzer <S in "P is the set of elements of P € "P 
such fliat PSP^f e S for aU S G S. Note fliat S C Af. 
Elements of the normalizer preserve the stabilizer subspace, 
Si{N\-ilj)) = NSj\ip) = N\^) for N € J\f. The group 
Af/S is a Pauli group on k qubits, i.e. Af/S acts as encoded 
operators on the stabilizer code. 

If an element E of the Pauli group anticommutes with 
at least one stabilizer element Sk ({E,Sk} ~ 0) then this 
element takes a state encoded into the stabilizer subspace into 
one which is outside of this subspace. This can be verified 
by noting that for e Hs, E\iP) = ESk\ip) = -SkE\ijj) 
which tells us that E flips the sign of the eigenvalue of Sk- 
This imphes that if we encode quantum information into 
our stabilizer subspace and then a Pauh operator E which 
anticommutes with at least one Sk acts on the system, then 
we can detect this event by measuring the eignenvalue of 
Sk- Elements of the Pauh group which anticommute with at 
least a single Sk are elements of the Pauh group which are 
not in M. 

This leads to the following characterization of the capabil- 
ities of a stabilizer code to correct errors. Let {Ei} be a set 
of Pauli elements such that eJEj ^ AT — <S for aU j, j. Then 
{Ei} is a set of correctable errors for this stabilizer code. 
The weight of the smallest E^Ej which is in — iS is the 
distance d of the stabilizer code. A quantum code which can 
correct up to t errors must have a distance at least 2t+l. A 
quantum error correcting code which uses n qubits, encodes 
k qubits, and has a distance d is called a [[n, k, d]] quantum 
code. 

II. Generalized Shor Codes 

Here we review codes which are simple generalizations 
of Shor's original [[9,1,3]] quantum error correcting code. 
The codes we construct later in this paper will be subsystem 
versions of these codes, achieving the same error correcting 
properties as a generalized Shor code, but with significantly 
reduction in the number of stabilizers which need to be 
measured in order to perform quantum error correction. 

To construct generalized Shor codes we first consider what 
a classical hnear code looks hke in the stabilizer formalism. 
Let C be a classical linear code with generator matrix G and 
parity check matrix P. Suppose that we wish to construct 
a quantum error correcting code which corrects bit flips 
(X errors) just as this classical code corrects bit flips. The 
stabilizer for this code is then obtained simply from the rows 
of the parity check matrix. In particular we can construct the 
StabiUzer generators 

n 

Si = (^Z''^'i. (1) 



Clearly these stabilizer elements all commute with each other 
because they are all made up of either I or Z Pauli operators, 
and they are independent of each other since the rows of 
Pi j are linearly independent. It is clear that measuring these 
stabilizer generators is equivalent to using the parity check 
matrix to calculate the syndrome for this stabihzer code. 

What are the logical operators for this code? The logical 
X operators for the k encoded qubits can be obtained directly 
from the generator matrix, G. In particular we can choose 
the ith logical X operator to be 

n 

X^^I^X^^-. (2) 

i=i 

These operators are certainly not in the stabilizer since they 
are made up entirely of X operators and they commute with 
all of the elements of the stabilizer because 

n n 
k=l 1 = 1 

^ {-l)^"^=iG.->'P^->'SjX,^S,X, (3) 

where we have used the fact that X^ILi Gi,kPj.k ~ 0. 
Further we see, as expected, there are k of these encoded 
X operators. What about the encoded Z operators? Well it 
is always possible to construct these operators from tensor 
products of Z and / operators. In particular we let the ith 
logical Z operator be 

n 

Z.^^Z^l^. (4) 
i=i 

What is -Pfj? Let S be the subspace of (F2)" spanned by 
the rows of the parity check matrix P. Then is a fc by n 
matrix whose row vectors are linearly independent from the 
rows of P, are themselves linearly independent, and which 
together with the rows of P span the entire space (F2)". 
In other words P'^ is made up of rows which together with 
the rows of P form a basis for the full space (F2)". The 
choice of a generator matrix G forces a particular choice 
(up to row multiplications) for P^ such that it is a logical 
Zi operator. Encoded Y operators can be obtained directly 
from Yi — iXiZi. 

In a similar vein to the construction of P"^, we can define 
a n — fc by n matrix G^ such that the rows of this matrix 
are linearly independent of the other rows of G, linearly 
independent of each other, and together with the rows of G 
form a basis for (F2)". We can then define the operators 

n 

E, = ^X'^l^. (5) 

Now the set of detectable bit flip errors is easily defined. It 
is any error which can be expressed as a product of at least 
one Ei and any number (including zero) of Xi operators. The 
reason these are detectable is that products of Ei operators 
are guaranteed to anticommute with at least one Si. 

Table m below presents a useful way of thinking about the 
structure of classical linear codes in the stabihzer formalism. 



In this table, each box shows the matrices used to construct 
a group for the linear classical code and the number of 
independent generators of the group is listed at the top of 
each column. In addition to the stabilizer and encoded X 
and Z operators, another group is formed from G^. We have 
labelled this latter group the group of "pure errors." These 
are errors which do not affect the encoded information and 
are detectable errors. 



TABLE I 

Classical linear codes in the stabilizer formalism. 





n — k Generators 


k Generators 


'I'cnsor product of 
Z and / operators 


P (Stabilizer) 


(Encoded Z) 


Tensor product of 
X and / operators 


(Pure EiTors) 


G (Encoded X) 



Finally let us note that the above construction can be 
used to correct phase flip errors instead of bit flip errors. 
This can be done by interchanging the Z Pauli operators 
in the stabilizer with the X Pauli operators. Replacing Z 
Pauli operators in the logical Z operator with the X PauU 
operators will produce a new encoded Z operator for this 
code. Replacing X Pauli operators in the logical X operator 
with the Z Pauli operators will produce a new encoded X 
operator for this code. 

We are now ready to describe generalized Shor codes. 
Shor's basic idea was that one could construct a quantum 
error correcting code for errors on qubits by concatenating 
a code designed to deal with bit flips (X errors) with one 
which is designed to deal with phase flips {Z errors). Notice 
that there is an asymmetry in this construction: one should 
decide which of these errors to be dealt with on the lowest 
level of the concatenation and which should be dealt with at 
the second level of the concatenation. 

Let Ci and C2 be two [ni,ki,di] and [n2,k2,d2] linear 
codes. Let Ci and C2 have generator matrices, Gi and G2 
respectively, and have parity check matrices Pi and P2 
respectively. In generalized Shor codes we use these two 
codes to construct a [nin2, fcifc2, min((ii, c?2)] quantum error 
correcting code. To do this we proceed as follows. Take nin2 
qubits and partition them into 712 blocks of size ni. For each 
of these blocks of ni qubits we can define a quantum error 
correcting code from Ci which is designed to correct bit flip 
errors as above. Each of the rii blocks encodes ki qubits. We 
can now use these encoded qubits to protect against phase 
flip errors. The first observation we need is than any tensor 
product of Z and / acting as an error on one of our blocks 
will act as a phase error on one or more of the encoded qubits 
times an element of the stabilizer. The reason for this is that 
every tensor product of Z and / is either in the stabilizer 
of quantum error correcting or is a product of a stabilizer 
operator and encoded Z operators. This allows us to consider 
phase errors on blocks as phase errors on the encoded qubit. 
To protect against phase flip errors we begin by picking the 
ith encoded qubit from every block. We can then use each of 
these from all blocks to construct a quantum error correcting 



code for phase flips using the C2 code. Since for each choice 
of i we obtain a code with k2 encoded qubits, this code can 
be used to store kik2 qubits. Further, the code will be of 
distance di for bit flip errors and ^2 for phase flip errors. 
For Y errors the code will be of distance niin(di, ^2). Thus 
the distance of the code will be min((/i, ^2). 

Generalized Shor codes are easy to construct, but do not 
have nice asymptotic error correcting properties. However, 
they are conceptually extremely easy to understand and 
their error recovery routines are easily obtained from their 
constituent classical linear code recovery routines. One im- 
portant point about these codes, however, is that they have 
an asymmetry in their recovery routines due to the choice 
of whether to concatenate bit flip codes with phase flip 
codes or vice versa. In particular for each the lowest level 
of concatenation, error correction should be performed for 
every single block. Indeed it is easy to see that the stabihzer 
for the generalized Shor code constructed above is generated 
by a set of (ni — fci)n2 + (n2 — ^2) independent operators. 
We wiU show below that using the notion of a subsystem, 
this can be reduced, for every generalized Shor code, to 
(ni - ki)k2 + {n2 - ^2)^1. 

III. Quantum Error Correcting Subsystem Codes 

In stabilizer codes one encodes information into a sub- 
space of a quantum system. However the most general way 
to encode information is not to encode it into a subspace, 
but instead to encode it into a subsystem[15]. Let us briefly 
review this concept and describe the notion of quantum error 
correcting subsystem codes. 

Suppose we have a Hilbert space Tl. Then one method for 
encoding quantum information is to encode this information 
into a subspace of H. In particular if H is the direct sum 
of two subspaces, Tic and Hd, "H = Tic © TId, then we 
can encode quantum information into one of the subspaces. 
He- In addition to the notion of a direct sum, ®, of two 
Hilbert spaces, another notion for combining two Hilbert 
spaces is to construct the tensor product of these two Hilbert 
spaces. Thus, for example, we can combine two Hilbert space 
Tic and Hd as 7Y = He <E) Hd- Then we can encode 
quantum information into the subsystem Tic- Notice that 
such an encoding, for a fixed encoding into Hn^sa subspace 
encoding, but, without such a specification, the encoding is 
not a subspace encoding. 

By repeatedly constructing subsystems and subspaces on 
Hilbert spaces, we can, most generally decompose a Hilbert 
space into a multiple direct sum of multiple tensor products 
of Hilbert spaces (since the process of direct sum and tensor 
product obey a distributive law.) Further if we single out a 
single one of these Hilbert spaces, call it He, then we may 
coUect the other terms in such a decomposition so that the 
Hilbert space decomposes as 

H = {Hc® Hd) e He. (6) 

The decomposition described above, which is the most gen- 
eral for encoding a single Hilbert space, can be described as 
taking a Hilbert space H and decomposed it into a subspace 



He and a perpendicular subspace, H^. On this perpendicular 
subspace we have further decomposed this into a tensor 
product of two subsystem Hilbert spaces, = He ®Hd- 
Thus if we are going to encode quantum information to the 
subsystem He we can do this by preparing the quantum state 

P = (pc <8> Pd) e Ob (7) 

where pc is the density matrix of the encoded quantum 
information, is information encoded into the subsystem V 
(which can be arbitrary) and 0^ is the all zero matrix on the 
subspace He- this point we can see one of the particular 
features of subsystem encodings: if we act nontrivially on 
the subsystem "D then the quantum information encoded into 
the subsystem C is not affected. In other words, information 
encoded into a subsystem is not affected by information 
encoded into different subsystems. Encoding into subsystems 
has been used most notably in noiseless subsystems[15], 
[26], [27] and communicating without a shared reference 
frame[28], as well as being essential to an important trans- 
form in quantum information theory, the quantum Schur 
traiisform[29], [30]. 

What does this mean for the theory of quantum error 
correction? Suppose that we encode quantum information 
into a subsystem He of a Hilbert space which has been 
decomposed as W = {Tie ® Hd) ® He- Next suppose 
that a quantum operation, corresponding to some quantum 
error, occurs on our system. Then the goal of quantum error 
correction is to restore the information encoded into the 
subsystem He- In the case where we have a subspace code, 
i.e. when Hd = ^, then we must apply an operation which 
correctly restores the quantum information encoded into the 
subspace He- If on the other hand we have a subsystem 
code, Hd 0, then we must apply an operation which 
correctly restores the information encoded into the subsystem 
He, but we do not care what happens in this procedure to 
the information encoded into Hd- other words, if we are 
to perform quantum error correction on a subsystem code, 
then the error recovery routine need only correct the error 
modulo the subsystem structure. We need not be worried 
if information encoded into Hd is destroyed by the entire 
error/recovery routine, as long as the information in He is 
correctly restored. 

Suppose that we wish to protect our quantum information 
from a set of errors Ea after we have encoded the quantum 
information into a subsystem as described above. Suppose 
that \i) ® |j) is a basis for the subspace He <8) Hd, then 
a necessary and sufficient condition[16], [19] for the set of 
errors Ea to be correctible is that 

<S> {k\)ElEf,{\j) |/» = SijCafi. (8) 

Notice that this condition does not depend on the \k) and 

10. 

For a more complete description of quantum error correct- 
ing subsystems we refer the reader to [16], [17] which details 
not just the notion of a quantum error correcting subsystem, 
but also the notion of a operator quantum error correction. 



which is a complete method for deaUng with quantum error 
correcting subsystems. 

IV. Subsystem Codes from Two Linear Codes 

We now turn to the construction of a new class of quantum 
error correcting subsystems. We will begin by detailing the 
construction and then proving that our construction has the 
error correcting properties which we claim. Our construction 
follows, in rough outline, that presented in [21]. 

A. Subsystem Code Construction 

Suppose we are given two classical linear codes, Ci and 
C2, which are a [ni,fci,di] code and a [712,^2,^2] code, 
respectively. We wiU now show how to use these codes 
to construct a quantum error correcting subsystem code 
which is a [[nin2, A:ifc2, inin((ii, ^2)]] code. Let Pi and P2 
denote the parity check matrices and Gi and G2 denote the 
generator matrices for the two codes Ci and C2 respectively. 
From the rows of Pi we can construct a stabilizer code on 
rii qubits. In particular we can define the rii — ki stabihzer 
operators Si = ^^^^^^'^^''^ Call the stabilizer group 
generated by this set of operators .Si = (^i, . . . , Sn-^-ki)- 
In a similar manner we can use P2 to construct a stabihzer 
code. In particular define the n2 — k2 stabihzer operators 
Ti = (8)"|,iX(^2)., This forms a stabihzer group ^2 = 

(Tl, . . . ,T„2_fe2). 

These codes are classical codes except that the second 
code is not in the computational basis |0), |1), but is instead 
in the dual |+), |— ) basis. The first code is designed to cor- 
rect |_ ''^~"'" J bit flip errors (Pauli X errors) while the second 
code is designed to correct [^^^^J phase flip errors (Pauli 
Z errors.) For both codes, we can follow our construction of 
classical codes and construct encoded operators. For the first 
code call the encoded X and Z operators (Xi)i and {Zi)i 
and for the second code call the encoded X and Z operators 
(X2). and {Z2)^. 

Now we will show how to produce a subsystem code using 
these codes. Put n = nin2 qubits on a rectangular m x n2 
lattice (this lattice is for illustrative purposes only, and is not 
a necessary part of the code.) We will now use the stabilizer 
codes operators Si in the columns and the stabilizer code 
operators S2 in the rows to construct a nonabelian group T 
on these n? qubits. More specifically, let 7i be the stabilizer 
group made up of letting Si operators acting on all of the 
columns of the lattice and let T2 be the stabihzer group made 
up of letting 1S2 operators acting on all of the rows of the 
lattice. Then the group we are considering, T, is the group 
generated by the elements of Ti and T2. 

The group T is clearly nonabelian. From T we can further 
construct an abelian invariant subgroup (invariant meaning 
that all elements of T commute with the elements of the 
subgroup.) To do this, we do the following. Take one of 
the stabihzer operators from <Si. Take one of the codewords 
from C2, call it v. Now construct an operator on our n qubits 
which has acting on each column j (where S*" = /.) 
Clearly these elements are in T. Further they commute with 
all of the elements of T since in a particular row they are 



made up of Z operators which have the form ^^l-^Z'"^ 
and in a particular column they are made up of elements 
of <Si. We can similarly take stabihzer operators form S2 
and codewords w from Ci and construct operators which 
act like acting on each row j. These operators will also 
commute with all of the elements of T. So now we can form 
our abehan invariant subgroup <S as the group generated by 
these stabilizer operators for aU possible stabilizer codeword 
combinations in both rows and colunms. 

So now we have a structure set up where we have a 
non-abelian group T and an abelian invariant subgroup of 
this group, S. There is another set of operators which are 
important, which wiU correspond to the logical operators on 
the code £. Suppose we take an encoded X operator for 
the stabihzer code Si, call it {Xi)i, and take an encoded 
operator for the stabilizer <S2, call it (^2)^. Then we can 
form an operator Xij acting on our qubits by putting 
{Xi)i in each column j where {X2)j acts non-trivially as 
X. Similarly from {Zi)i and {Z2)j we can construct an 
operator {Zij) which is {Zi)i in each column j where {Z2)j 
acts non-trivially as Z. It is easy to see that Xij and Zij 
commute with the group T since in each row or column 
they look hke encoded operators. Further these operators 
also anticommute with each other if and only if their indices 
match, {XijjZk^i} = 6i^kSj,i, and hence commute with each 
other otherwise. These operators further form a group which 
is isomorphic to a Pauli group on kik2 qubits. 

Next we need to discuss how to put T, S and C together 
to form a subsystem code. This is nearly identical to the 
procedure described in [21]. 

To do this, it is convenient to adopt explicit forms for the 
operators we have described above in a simple notation. Let 
M be a ;/ by n matrix with entries either or 1. Then we 
define as the operator on our n by n qubits which is a 
tensor product of P and / operators which acts on the qubit 
at the ith row and jth column as P*^» 3 . In this notion it is 
easy to see that every element of the Pauli group on our 
qubits can be expressed as i'^Z^X^ where k € {0, 1,2,3} 
and A and P are n by n 0/1 matrices. We will now proceed 
to use this notation to express the operators in T, S and C. 

First consider T. Consider operators in T which have 
elements of Si which lie in a column and are products of 
Z operators. In our new notation, these operators can be 
written as Z^ where Aij = {Pi Pi)iSjjQ where pj is a 
ni — ki binary row vector and jo is the column where this 
operator acts. Similarly, we can construct operators from ^2 
lie in a row and are products of X operators. They can be 
expressed as X^ where Bij = Sis„{p2 P2)j where io is the 
row where this operator acts and is a 71,2 — fc2 binary row 
vector. In order to make a distinction which will be useful 
later, it is useful to express the delta functions in the above 
expression as follows. Since the rows of Gi and form 
a basis for the entire space (F2)" we can express Sjjg as 
{gfGi + {glYGDj for some choice of length k binary row 
vector and length n — k binary row vector (.gj)"^. We 
can perform a similar decomposition for the X operators. 
Since multiphcation now corresponds, up to a phase factor, 



to addition, it is then easy to see that every element of the 
group T can be expressed as t{Q, Q'^, R, R^,p) = i^Z^X^ 
with 

A = P[{Q)G2 + P^iQ'')G^^ (9) 

and 

B = Gl{R)P2 + {G'if{R')P2 (10) 

where Q is a ni — fci by fc2 0/1 matrix, is a ni — ki 
by n2 — A;2 0/1 matrix, R is a ki by n2 — A;2 0/1 matrix, 
R'^ is a ni — ki by n2 — ^2 0/1 matrix, and p e {0, 1, 2, 3}. 
Note that we have extended the group slightly by adding 
in a phase factor of i. This is for convenience sake when 
describing certain Pauli subgroups of T. 

Having described T in our notation, we now turn to S. 
5 is a subgroup of T and thus we can express it again 
as a t{Q,Q^,R,Rc,p). In particular from the definition of 
the S we see that its elements are of the form s{Q,R) = 
t{Q, 0, R, 0, 0), i.e. they are elements of T with = 0, 
i?^ = 0, and 39 = 0. 

Next we express elements of C in this notation. Elements 
of £ are formed from encoded operations from the two codes. 
In particular it is easy to see that they can be expressed as 
l{U,V,p) =iPZ^X^ with 

A = {P^f{U)G2 (11) 

and 

B = G{{V)P^ (12) 

where U and V are ki by k2 0/1 matrices and p € 
{0, 1, 2, 3}. Again we have added an extra phase factor of i. 

Finally let us notice that every element of the 
Pauh group on our -n? qubits can be expressed as 
o{Q, Q", R, R", U, U", V, V'=,p) = iPZ^X^ with 

A = Pl{Q)G2+PTmGl+{Pif{U)G2+{Pi)'^{U'')Gl 

(13) 

and 

B = Gl{R)P2 + {G'if{R'=)P2+Gl{V)P^+{Glf{V'')P^ 

(14) 

where Q, Q'^, R, i?^, U, and V have the dimensions listed 
above, is a ki by n2 — 0/1 matrix, V is a ni — ki 
by k2 0/1 matrix, and p e {0, 1, 2, 3}. 

We can now define our subsystem code. We begin with 
<S. <S is an abehan subgroup of the Pauli group which does 
not contain — . Thus it is a stabilizer group. From this 
stabilizer group we can from a stabilizer code S. Operators 
which are in T and £ are then in the normahzer of this 
stabiUzer code, since all of the elements of T and C commute 
with all of the elements of S. Further we can think of the 
operators from C and T acting on different logical qubits for 
the stabilizer code, since all of the elements of C commute 
with all of the elements of T. 

The subsystem code can now be defined. Since <S is a sta- 
biUzer code, we can label subspaces of the qubits Hilbert 
space by the ±1 eigenvalues of the stabilizer generators. 
There are (ni — k\)k2 such stabihzer generators made up of 
tensor products of I and Z operators and fci(n2 — ^2) such 



StabiUzer generators made up of tensor products of / and X 
operators. Thus the number of generators for this stabiUzer 
group is (rii — k{)k2 + k\ (n2 — fe)- Next notice that C forms 
an encoded Pauli group acting on kxk2 qubits. To see this 
notice that l{U,v,p) follows the multipUcation rules of a 
Pauli group as if for kik2 independent qubits. Now consider 
T/S, the group T after we divide out the stabilizer group. 
This group is made up of operators t(0, Q'^, 0, This 
group is an encoded Pauli group on (m — fci ) (n2 — A;2 ) qubits. 
Putting this together we have a stabilizer with (ni — ki)k2 + 
ki{n2 — ^2) generators, encoded T/S operators which act 
as a Pauli group on (rii — fci)(n2 — ^2) qubits, and encoded 
C operators which act as a Pauli group on kik2 qubits. The 
total of these generators and the number of encoded qubits 
is {rii - fci)(n2 - ^2) + {rii - ki)k2 + ki{n2 - fc2) + fcifc2 = 
nin2- This implies that the T/S encoded operators and 
the C encoded operators form an exhaustive Ust of encoded 
operators for the stabiUzer code <S. 

So how do we define our subsystem code? Since there are 
(ni ~ki)k2 + ki (n2 — ^2) stabilizer generators, we can label 
subspaces of dimension 2"i"2-(ni-fci)fe2+fei(n2-fc2) 

±1 eigenvalues, call them Si, of these stabilizer generators. 
Call s the (ni — ki)k2 + ki{n2 — ^2) tuple of these values. 
Further for each such subspace there is now a tensor product 
between encoded logical operators from C and those from 
T/S. Thus we can find a basis for our Hilbert space on nin2 
qubits such that it decomposes as 

n= wf/^®Hf (15) 

se{0,l}<"l-'=l)'=2+'=l('>2-*:2) 

where dimWf = 2''^''^ and dirnHj = 2("i-'=i)("2-'=2). It 
is now useful to describe how elements of S, T/S and £ 
operate on this decomposition. Elements of S act as either 
±1 on each subspace. In particular they act as 

e, , , ('>l-'=l)'=2 + ''l("2-''2) „. ^ _ ^ 
(-l)^i=l ^'I^I 

S6{0,1}("1 ^'"l'''2+''l("2-'=2) 

(16) 

Elements of £ act as encoded qubits on the Hf subsystems 
/0L(s) (17) 

while elements of T act as encoded qubits on the Hj 
subsystems 

T(s)0 7. (18) 

Our subsystem code can now be defined. We will encode our 
quantum information into the all s, = +1 subspace and the 
corresponding subsystem. This encoding will encoded 
kik2 qubits and the logical PauU operators on this code come 
from £. Note that elements of T can always be expressed 
as operators which do not act on this subsystem. This is 
a subsystem degree of freedom which makes our code a 
quantum error correcting subsystem code. 



B. Subsystem Error Correcting Routine 

Having identified the subsystem we are encoding into and 
the representation theoretic structure of operators in S, T 
and £ we can now turn to the error correcting procedure for 
this code. We will show that this code can detect single qubit 
errors of weight mm{di, c?2) and describe the error recovery 
routine for this code. This error recovery routine corrects the 
information in our subsystem but may act nontrivially on the 
subsystem Hj . 

Suppose that Ci and C2 can correct the sets of errors £1 
and £2 respectively. We will now show how this allows us to 
correct bit flip and phase flip errors modulo the subsystem 
structure. 

First consider an X error on our code. Suppose that ei is 
a correctable error for code Ci. Let Ei be a ni by 712 matrix. 
We will show that every error with Bij = {ei)i{Ei)ij 
is a correctable error for our subsystem code. To see this we 
first note that if we express this error as in Eq. il4l as 

B = G'[iR)P2 + iG'lf{R')P2 + G{{V)P^ + {G'lfiV^)P? 

(19) 

then we can turn this into a product of an element of T and 
one which is not in T, X^^X^^ where 



Bi = Gi {V)P^ + (G^)' (V")P; 



2 ' 



and 



(20) 



(21) 



Since X^^ is an element of T it cannot act as an error on 
the information encoded into our subsystem. Therefore we 
can consider the error to be purely of the form X^^ . Now 
since the errors we are considering have B matrices of the 
form X^ with Bi j = {ei)i{Ei)ij this means that we can 
restrict the rows of Ei to be from the subspace spanned by 

pc 

Now suppose that we measure the elements of S made up 
of tensor product of / and Z operators. These operators are 
of the form Z"^ where 



A = iPifiU)G2 



(22) 



P| and G2 act as encoded X and Z operators, respectively, 
for the quantum version of the code €2- Thus by measuring 
the stabilizer generators which are tensor products of / and Z 
operators, we can, for each encoded qubit in the second code, 
make a measurement of the Pi for these encoded qubits. If 
the error ei is a correctable error for Gi, then we can apply 
the appropriate eipj operators where pi are the appropriate 
row vectors from P2 . The effect of this correction procedure 
will be to restore the information encoded into the subsystem 
up to the operator X^'^ which is an element of T. Thus we 
see that, as claimed, that we can correct errors which are of 
the form X^ with Bij = {ei)i{Ei)ij. 

A similar prescription applies for Z errors by measuring 
the stabilizer generators which are made up of tensor prod- 
ucts of / and X operators. If 62 is a correctable error for the 
code C2, then this will correct errors of the form Z^ with 
A-ij = (62) j (£'2)1 J - If these two procedures are carried out 



one after another they will also correct errors which Z"^X^ 
with A and B correctable as above. 

Thus we see that by measuring the generators of S we 
can correct errors related to the original codes Ci and C2. 
What is the distance of this code? For bit flip errors the 
distance will be di since we the smallest error of the form 
X^ with Bij = {ei)i{Ei)ij has only single X errors in a 
row. Similarly the distance for phase flip errors will be c?2- 
The full distance must include Z^X^ errors and thus the 
distance is min((ii, ^2)- 

C. Savings Over Generalized Shor Codes 

In a generalized [[nin2, fcifc2, inin((ii, ^2)]] Shor code we 
have seen that error correction is achieved by measuring 
(ni — fci)n2 + (n2 — fc2) stabilizer generators. In our construc- 
tion of a [[nin2, fcifc2, min((ii, ^2)]] subsystem code above 
we have achieved the same parameters for the code but now 
using only (ni — fci)fc2+fci("-2 — ^2) stabilizer measurements. 
Both methods use two classical linear error correcting codes 
to construct a new quantum error correcting code. Indeed, 
the subsystem codes we present are nothing more than 
generalized Shor codes with certain stabilizers which do 
not add to the error correcting distance removed[20]. The 
subsystem code versions of generalized Shor codes have 
considerable advantages over the subspace code when it 
comes to the complexity of the error recovery routine, 
providing a quadratic savings in the number of stabilizers 
which need to be measured. 

V. Examples 

In this section we present a few examples of our code 
construction. 

A. Redundancy Code 

This is the construction presented in [21]. Let Ci and C2 
both be a simple n qubit redundancy code with generator 
matrix G = (1,1,....!) and parity check matrix 



P = 



1 1 
1 

••• 



... 
... 

1 1 



(23) 



This redundancy code is a [n, code. The resulting 
subsystem code is a [[n^, 1, n]] code. Properties of this code 
are described further in [21]. Notice that a generalized Shor 
code constructed from this redundancy code requires the 
measurement of rt^ — 1 stabilizer operators. Our subsystem 
code achieves the same parameters for the code but using 
measurements of only 2{n — 1) stabilizer operators. This 
increase in efficiency combined with other nice properties of 
this code has recently been shown to improve the threshold 
for fault-tolerant quantum computation[24]. 



B. A [[49,1,5]] Code Which Outperforms a Concatenated 

Steane Code 

Consider using the Hamming [7, 4, 3] code for the codes 
C\ and C2. Our construction will yield a quantum error 
correcting subsystem code with parameters [[49, 16, 3]]. Now 
consider using those 16 encoded qubits in a redundancy 
[[16,1,4]] code as described in the last subsection. If we 
use the error recovery for our [[49, 16, 3]] code followed by 
error correction for the [[16,1,4]] code this will produce 
a code which can correct arbitrary 2 qubit errors and is 
thus effectively a [[49, 1, 5]] code. The number of stabilizers 
which need to be measured for this use of these codes 
is 24 + 6 = 30. Another option is to use on 9 of the 
16 encoded operators in the redundancy subsystem code 
described above (allowing any error to occur on the other 
subsystem.) In this case one achieves a subsystem code with 
parameters [[49, 1, 5]] but with a number of stabilizers given 
by 24 + 4 = 28. 

This should be compared with the normal (not optimal) 
use of concatenating the Steane code [[7,1,3]] code[31] 
with itself. There one uses error correction on the Steane 
code [[7, 1, 3]] for multiple levels of the concatenation. This 
results in a code which is effectively a [[49,1,5]] code. 
This concatenation scheme will require the measurement of 
measure 42 + 6 = 48 stabilizer operators. Thus we see 
that a considerable savings of 18 (or 20) less stabiUzer 
measurements. 

VI. Conclusion 

We have shown that for every generalized Shor code 
there is an subsystem code with the same parameters but 
which requires significantly fewer stabilizer measurements 
in order to perform quantum error correction. These codes 
are generalization of the codes presented in [21J and are 
in the class of stabilizer subsystem codes described in [20]. 
Recently Aliferis and Cross [24] have used the subsystem 
codes described in [21] to significantly improve the provable 
threshold for fault-tolerant quantum computation. A major 
open question is whether the subsystem codes described in 
the work described here can lead to similar and perhaps 
greater increases in the threshold for fault-tolerant quantum 
computation. 
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